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ABSTRACT 

We present broad band analysis of the Be/X-ray transient 1 A 1 118-61 by Suzaku at the peak of its 3rd ob- 
served outburst in January 2009 and 2 weeks later when the source flux had decayed by an order of magnitude. 
The continuum was modeled with a cutoff pi model as well as a compTT model, with both cases requiring 
an additional black body component at lower energies. We confirm the detection of a cyclotron line at ~ 55 ke V 
and discuss the possibility of a first harmonic at ~ HOkeV. Pulse profile comparisons show a change in the 
profile structure at lower energies, an indication for possible changes in the accretion geometry. Phase resolved 
spectroscopy in the outburst data show a change in the continuum throughout the pulse period. The decrease in 
the CRSF centroid energy also indicates that the viewing angle on the accretion column is changing throughout 
the pulse period. 

Subject headings: X-rays: stars — X-rays: binaries — stars: pulsars: individual (1 A 1 118-61) — stars: mag- 
netic fields 



1. INTRODUCTION 

The Be/X-ray binary transient 1A 11 18-61 was serendip- 
itously discovered during an observation of the nearby bi- 
nary system CenX-3, when an outburst was detected in De - 
cember of 1974 by the Ariel-5 satellite dEvles et all 1 19751) . 
A second, similar outburst occured in January of 1992 and 
was observed by the Burst and Transient Source Experiment 
on the Compto n Gamma Ray observatory CGRO/BATSE 
(ICoe et al.l [i994). The measured peak flux was ~ 150 mCrab 
for the 20-100keV energy range, similar to the 1974 out- 
burst. The source showed an elevated emission throughout 
the next ~ 30 d ays after the decay of the main outburst (see 
ICoe et all 1 1994 Fig. 1). The third and most recent outburst 
occurred on 2009, January 4 and was detected by the Swift 
Burst Alert Telescope BAT (Man gano etal.ll2009t iMangano 
120091) . It reached a peak flux of ~ 500 mCrab in the 15-50keV 
energy band. This last outburst was monitored with Swift and 
the Rossi X-ray Timing Explorer (RXTE) as well as with two 
long Suzaku pointings and one observation with INTEGRAL 
during a flaring episode ~ 50 days after the peak of the main 
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outburst dLevder. Walter & Lu binski 2009|). 

Pulsations with a period of 405.3 ± 0.6s were ob- 
served during the 1974 outburst and were initially at- 
tributed to the orbital period of two compact objects 
dives. Sanford & Bell Burnetii 119751) . Shortly afterwards it 
was suggested that the period stems from a slow rotation 
of the neutron star (NS) itself (Fabianl [l975l) . During the 
1992 outburst pulsations with a period of ~ 406.5 s were de- 
tected up to 100 keV, showing a broad, asymmetric, single 
peak pulse profile above the lowest BATSE energy of 20keV. 
The pulse period decreased throughout the decline of the out- 
burst with a rate of ~ -0.016 s/day and it appeared constant at 
~ 406.35 s for the time of the elevated emission. During the 
2009 outburst a similar period evolution was observed with 
RXTE resulting in a pulse period of P sv i B = 407.719 s, and 
Pspm = -4-6 x 10" 7 s/s rj -0.04 s/day dDoroshenko et al.l2010l) . 
Furthermore, the lower energies showed a more complex 
pulse profile with two peaks below an energy of ~ lOkeV. 
Due to the short duration of the Suzaku observations with re- 
spect to the pulse period, the derived RXTE values are used 
for determining the pulse profile and phase resolved spectra 
in this paper. 

The optical counterpart was identified as the Be- star 
Hen 3-640/Wray 793 by iChevalier & Ilovaiskvl (fl975l) 
and classified as an 09.5IV-Ve star with strong Balmer 
emission lines indicating an ex t ended envelope by 
Uanot-P acheco. Ilova iskv & Chevalier! (1198 11) . The over- 
all spectrum was found to be similar to other known 
Be/X-ray transients, such a s X-P er and A 0535+26 
( Villa d~Giovannelli & Polcarol Q992, and references 
therein). The distance was estimated to be 5 ± 2 kpc 
dJanot-Pacheco. Ilovaiskv & Chevalier! IT981) and was con- 
firmed by ICoe & Paynd d!985l) . along with the spectral type 
classification, using UV observations of the source. EX- 
OS AT obs erved X-ray emiss ion from 1A 1118-61 between 
outbursts (Motch et al. 1988), thus indi c ating a continuous 
low level of accretion. iRutledge et al.l ( 120071) reported on 
pulsations in the low luminosity state observed with Chandra, 
making it only the third known HMXB transient after 
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A 0535+26 and 4U 1145-619 for which this behavior has 
been observed. 

A study of the H n emission line before and during the 
1994 outburst dCoe et all 1994) showed a strong correlation 
between its strength and the observed X-ray activity, indicat- 
ing the existence of a very large disk around the star. The 
fluctuation of the equivalent width of H a indicated possible 
instabilities in the disk, which were enhanced when the star 
passes through the periastrion. The analysis of the UV con- 
tinuum and line spectra has indicated that the photospheric 
emission from the Be star was not affected b y the X-ray radi - 
ation, similar to the case of A 0535+26 (Ide Looreet al.C l984). 

Until recently, the orbital period of 1A 11 18-61 was not 
measured and assumed values were of the order of 350 days, 
based on the P spin - P or \, relation (ICorbetl 119861) . and ~ 
585 days based on the equ i valent width of the H a line 
dReig. Fabregat & Cod fl997). Stau bert et all (120101) analyzed 
the pulse arrival time in RXTE data throughout the 2009 out- 
burst and established an orbital period of ~ 24 days with a 
very circular orbit around the Be-companion. 

Spectral fitting during the 1974 outburst indicated a vari- 
able power law index, where the lowest value was ~ 0.9 
duri ng the peak of the outburst, and ~ 1 .1 before a nd af- 
ter (lives. Sanford & Bell BurneU| [T975T) . iCoe et all (fl99l 
described the 1992 CGRO/BATSE data with a single- 
temperature, optically-thin, thermal bremsstrahlung (OTTB) 
with a temperature of 15keV. Due to the different energy 
ranges of the instruments and the fact that pulse-off-pulse 
data have been used, these results can not be directly com- 
pared with the 1974 o bservation. 

In the 2009 outburst Dorosh enko et al.l (120101) discovered a 
cyclotron resonance scattering feature (CRSF) at ~ 55 keV. 
CRSFs can be used to deduce the strength of the magnetic 
field at the pole region of the NS and have been observed in 
multiple sources at energies from ~ 15keV up to ~ 50keV 
(ICoburn et al.ll200l . 

Suzaku observed 1A 11 18-61 twice in 2009, once during 
the peak of the outburst and again ~ 13 days later when the 
flux had returned to a level slightly above the quiescent state. 
In this paper we will present the first detailed 0.7-200keV 
broad band spectrum for this source, confirming the detec- 
tion of a CRSF at ~ 55keV. The paper is divided as fol- 
lows: Section 2 describes the data accumulation and extrac- 
tion; Section 3 presents the broad band phase averaged spectra 
for the outburst (0.7-200 keV) and for the second observation 
(0.7-70keV), and discusses the CRSF and Fe-Line complex; 
Section 4 compares energy resolved pulse profiles for both 
observations and discusses the phase resolved analysis for the 
outburst data; finally Sections 5 and 6 present a discussion of 
our observations and summarize our findings. 

2. OBSERVATION AND DATA REDUCTION 

A sudden increase in activity of 1A 11 18-61 was de- 
tected with the Swif t IB AT instrument on January 4th, 2009 
(Mang ano et al.l I2009I) . The count rate increased steadily 
until January 15th when it reached the maximum value of 
~ 500mCrab and then decayed until January 27th where 
it showed a low emission leve l with periods of flaring 
dLevder. Walter & Lu binski 2009) until mid March and then 
returned to quiescence (see Fig. [TJ. Suzaku observed 
1 A 1 118-61 during the peak of the outburst on January 15th, 
2009 (MJD 54846.5, ObsID 403049010) with both of its 
main instruments: the X-ray Imaging Spectrometer (XIS; 
iMitsuda et al. 2007) and the Hard X-ray Detector (HXD; 
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FIG. 1. — Swift/BAT lightcurve during the 1A 1118-61 outburst, 
dashed lines mark the Suzaku observation as well as 2009 Jan 1. 
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Takahashi et al. 2007). A second observation was performed 
on January 28th, 2009 (MJD 54859.2, ObsID 403050010), 
~ 13 days after the main outburst, at the beginning of the pe- 
riod of elevated emission. 

Both observations were performed using the HXD nomi- 
nal pointing to minimize the pile-up fraction in the XIS in- 
struments and to enhance the HXD sensitivity for a possible 
CRSF detection. Of the four original XIS instruments only 
XIS 0, 1 and 3 were functional during the observing time. 
XIS and 3 are both front illuminated (FI) CCDs and XIS 1 
is a back illuminated (BI) CCD. To reduce pile-up, the XIS 
instruments were operated with the 1/4 window option with a 
readout time of 2 s and the burst option with only 1 s accumu- 
lation time for each readout cycle, reducing the exposure time 
for the XIS instruments by 50%. The data were taken in the 
3x3 and 5x5 editing modes which were then combined for 
the final spectral analysis. 

For the extraction the Suzaku FTOOLS version 16 (part 
of HEASOFT 6.9) was used. The unfiltered XIS data were 
reprocessed with caldb20090402 and then screened with the 
standard selection criteria as described in the ABC guideFT. 
Each detector and editing mode combination was extracted 
independently and individual response matrices and effective 
area files were created. For the final spectra the data from 
both FI detectors were combined (XIS and 3) and the re- 
sponse matrices and effective areas were weighted according 
to the accumulated exposure time of the different modes. The 
XIS data were grouped so that the minimum number of chan- 
nels per energy bin corresponded to at least the half width half 
maximum of the spectral resolution, i.e. grouped by 8, 12, 14, 
16, 18, 20, 22 channels starting at 0.5, 1, 2, 3, 4, 5, 6, and 
7keV, respectively (Nowak, priv. com). 

For the HXD data the Suzaku team provided the tuned PIN 
non X-ray background^ (NXB). Following the ABC Guide 
the cosmic X-ray background (CXB) was simulated and the 
exposure time was adjusted in both backgrounds by the pre- 
scribed factor of 10. The PIN data were grouped so that at 
least 100 events were detected in each spectral bin. 

GSO data were extracted using the FTOOL hxdpi with 
the newest gain calibration file from April, 2010. The NXB 
background files version 2.4 created by the Suzaku HXD in- 
strument team were used. The data were then binned to 64 
bins according to the Suzaku homepag4__ GSO data in the 

1 1 http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/ 

12 http://heasarc.nasa.gov/docs/suzaku/analysis/pinbgd.html 

13 http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/gso_newgain.html 
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range 70-200 keV were used in the spectral analysis. 

These selection criteria resulted in ~ 25 ks exposure time 
for each XIS instrument and ~ 49 ks for the HXD instruments 
in the first observation. The second observation had an expo- 
sure of ~ 2 1 ks and ~ 29 ks for the XIS and HXD instruments, 
respectively. 

2.1. Pileup correction 

For bright sources, such as 1A 11 18-61 
during the outburst, a strong pileup is ex- 
pected. Following the description provided at 
http://space.mit.edu/ASC/software/suzaku/ 
the S-lang routine aeattcor . si was used to improve the 
attitude correction file and the point spread function of the 
events. Then the tool pile_estimate . si was applied 
to produce a 2 dimensional map of the pileup fraction. The 
maximum values for pileup fractions were 10% and 15% 
for XIS0, 15% and 16% for XIS1, and 18% and 21% for 
XIS3 for the 3x3 and 5x5 editing modes, respectively. 
Regions with a pileup fraction above 5% were excluded 
from the extraction for each individual source event file. 
For the second observation the calculated maximum pileup 
fractions were < 5% and no regions had to be excluded for 
the extraction. 

3. PHASE AVERAGED SPECTRUM 

For the outburst observation we extracted broad band XIS 
spectra for 0.7-12keV (0.7-10keV for the BI XIS1), PIN 
spectra for 12-70keV and GSO spectra for 70-200 keV. For 
the second observation the GSO spectrum was not well con- 
strained and therefore was not included in the analysis. The 
final model included the Galactic and intrinsic absorption, a 
continuum that steepened at higher energies, an iron line com- 
plex, and a CRSF Excessive low energy flux was modeled 
with a black body component. In addition a lOkeV absorp- 
tion feature was required in the outburst observation. 

The Galactic column density was modeled with a single 
photon absorption (phabs) component, where the column 
density was confined between 1.1 and 1.4 xl0 22 cm -2 . The 
value determined by the NASA N H TooH for 1A 11 18-61 
is 1.22 x 10 22 cm" 2 . The A^h value was left free in the fits 
for the first observation. For the second observation the val- 
ues showed larger error bars with the cutoff pi model and 
had to be fixed when using a compTT model. The intrin- 
sic column density was modeled with the partially covered 
photon absorption model pcfabs to take the flux at lower 
energies (< 1 keV) into accoun t. All modeling was performed 
using the wilm abundances (Wilms, Allen & McCray 120001) 
with the vern cross-sections dVerner et al.lll996i) . The con- 
tinuum was modeled using a power law with an exponen- 
tial cutoff (cutof fpl). Using a power law with a Fermi- 
Dirac (FD) cutoff, one of the empirical models often applied 
to accretion powered pulsars, resulted in a cutoff energy of 
~ 0, making the FD cutoff effectively a cutof fpl, where 
the variable E cut actually reflects the folding energy of the 
model. The additional lOkeV feature was modeled using a 
negative, broad (a = 1 -2keV) Gaussian component. Such 
a feature has been previously observed in different sources 
and it is believed to ste m from an imprope r modeling of the 
continuum (see §6.4 in Cobu rnet al.ll2.002b . Mo re recent ex- 
amples of such a feature are, e.g. 4U 1907+09 dRivers et all 

14 http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 
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FIG. 2. — Broad band phase averaged spectrum of the outburst observation 
using the cutof fpl model. The residuals are shown without CRSF (b), 
with one CRSF at 55keV (c) and with the second CRSF at 112keV. The 
Inset shows the Fe line region and the best fit model for XIS1 and XIS 03. 

12010b and Cen X-3 dSuchv et all 12008b . Although this fea- 
ture sometimes is described as a broad emission line, e.g. 
ISuchy et"aTI (120081) . in this case a negative Gaussian line at 
^10keV fits best. When including this line in the overall best 
fit, the x 2 /dofs (degrees of freedom) decreased from 1173 / 
477 to the best fit value of 678 / 474, respectively. Due to 
the small differences in the instrument response, all detectors 
were coupled with a normalization constant, which was set to 
1 for the combined XIS 03 detectors. The final model is of 
the form const*phabs*pcf abs* (cutof fpl*gabs+ 
2*Gauss Fe ) +Gaussiokev+3*Gausscai- The best fit pa- 
rameters are mentioned in the text and partially summarized 
in TableQ] 

A second appr o ach fo r the continuum was initiated by 
iDoroshenko et a l. (2010) using the Compton ization model 
compTT developed by lTitarchuk & Hua ( 1995) instead of the 
cut of fpl model. The parameters of the compTT were 
again left independent between the FI XIS03 and the BI XIS 1 , 
but only small differences could be observed. The negative 
Gaussian at lOkeV was not necessary for this model. When 
comparing the spectral parameters in the outburst between 
these two models, i.e. CRSF, Nu and the spectral lines, to the 
previously described cutof fpl, the values were consistent. 

A comparison of the compTT values, such as the photon 
seed temperature To, the plasma Temperature kT e or the opti- 
cal depth t p showed consistency between the spectra of the 
outbu rst as seen by Suzaku and RXTE (IDoroshenko et al.l 
12010b . 
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table l 

Phase averaged spectral parameter. Model described in the text. 



Parameter 

phabs/V H [10 22 /cm 2 ] 
pcfabs /V H [10 22 /cm 2 ] 
covering fract. 
blackbody kT [keV] 
blackbody norm [10~ 3 ] 
Bout [keV] 

r 

A p i[10- 2 ]t 

^crsf [keV] 
ffCRSF [keV] 

TCRSF 

£FeK Q [keV] 

E FeKf3 [keV] 

Eq. Width K Q / [eV] 

Flux2-iokeV [10 36 ergs/sec ]* 

Cxisi/Cpin/Cgso 

X 2 /dofs 



Outburst" 

I 34+O.OI 
1,J -0.01 

I I 2+ 02 
66 +001 

u - -0.003 
7 ^+0.4 

'■ J -q,2 

35(1+0.022 
u.jju_ 0004 

c+o.n: 



'•SSL. 

58.2+™ 
14.H; 
60.1+5 9 

fi 41 4+O.OO6 

7 ia+0.03 
' ,1J -O.03 

51^/8+1 
8.79 

0.98/1.10/0.82 
678 / 474 



Cutoffpl 
2nd Obs. 
1.34 frozen 
1 3 5+ - 4 
071 +o.oi 

15+™ 

7 7+O.2 

-CUB 

3+3.0 

7+0 2 ? 4 
-0 04 
a+0.3 
8 -0.3 



Parameter 



-0..01 
26.3+ 

1.02 + 

5.8+ - 3 



58.2 frozen 
14.1 frozen 
15.4^ 

7.13 frozen 

1.72 
1.02/1.01/- 
574/431 



2nd Obs, free CRSF 
1 43+0.14 

13 6+ a5 

071 +o.m 

15+™ 

7 os+1.67 

24 7+'^ 

-J.s 

1 - ul -0.03 
5 7+ ' 3 



47.4+" 

5.7lf 

6.0+ 3 ? 

6 42+0 0' 
°-^ z -o.oi 

7.13 frozen 
S+ 2 
1.72 
1.02/1.01/- 
567 / 428 



46+!°/ 



compTT To [keV] 
compTT kT [keV] 
compTT r p 
compTT normt 



Outburst 

1 70+0.03 
1 - 2S -0.02 
12.00*^ 
,+0.01 

12.8+" 
4+0.04 

M 5 



CompTT 



0.71 



1.34 
7.6 



6.00+ 



n+o.w 

-0.18 
n n fi +0.003 

54 - 5 -U 

10 3+ 3 '° 
-0.1 

23 6+ 10 ' 4 
ZJ - -6.9 

6 42+0- 01 
°- -0.01 

7 13+0.03 
'■ -0.03 

51^/10! 3 
9.15 
0.99/ 1.18/1.10 
680 / 474 



2nd Obs. 
1.28 frozen 
15 1+ 01 
76+ 001 

it"" 51 

u -"+o.oi 
2 89+ 003 
44+0.06 

182+ff 

0.68+1; 33 

002+° 004 
u.uuz_ 001 

54.5 frozen 

10.3 frozen 

27 1+ 3 ' 2 

6 42+ 01 
°-^ z -0.01 

7.13 frozen 

4?-i8/ii-2 

1.72 
0.95/ 1.00/- 
567 / 428 



NOTE. — " Best fit values including first harmonic at ~ 1 10 keV, see text for details, t Units in Photons keV cm " s , | unabsorbed flux using a distance of 5 kpc 



10° 



> 

lo- 




ur 



10" 



nr 



10" 5 
6 
4 
2 

X 

-2 
-4 
-6 
6 
4 
2 

-2 
-4 
-6 



1.5 



1.0 



0.5 





\r XIS 03 




XIS 1^* 




6.0 6.2 6.4 6.6 6.S 7.0 7.2 



X 2 605/ 429 dofs 



X 557 V 426 dofs 



Energy [keV] 



FIG. 3. — Broad band phase averaged spectrum of the second observation 
using the cutoffpl model. Only XIS and PIN data were used. Residuals 
without a CRSF (b) and with one CRSF at ~ 47 keV. The Inset shows the Fe 
line region and the best fit model for XIS1 and XIS 03. 

When comparing the outburst spectrum with the second ob- 
servation, significant changes in the continuum parameters are 
observed. In the cutoffpl model the power law index soft- 
ened from 0.350!<j$g to 1 .01_g;g|. The cutoff energy £ cut in- 
creased from 18.7^| keV during the outburst to 24.7_| gkeV 
in the second observation. The temperature of the additional 
black body decreased from 0.51 +0 ' : 



, . 003 keV to 0.15!°;»lkeV 



after the outburst. The intrinsic column density and cov- 

+0.2 x 10 22 C m- 2 and 

22 „„,-2 



ering fraction increased both from 1 1 



<,+0.5 
-0.5 



1+0.01 
-0.01 



£,+0.01 

J -0.01 

in the second observation 

In the compTT model the changes between the two obser- 
vations were similar. The intrinsic column density and black 
body values were consistent with the values obtained with 
the cutoffpl model, where the intrinsic column density is 
12.0_g;_ x 10 22 cm" 2 and 15.1_^ x 10 22 cm" 2 and the cover- 



-0.7 



ing fraction is 0.71_§;§| and 0.76+o oj for the outburst and the 
second observation, respectively. When establishing the er- 
rors for the other spectral parameters of the continuum, these 
values had to be frozen to avoid a drift into an unreasonable 
part of the parameter space. The black body temperature de- 
clined from 0.53^0 {J} keV to 0.17_o(j] keV, similar values as in 
the cutoffpl model. The compTT values changed signif- 
icantly between the two observations. The photon seed tem- 
perature decreased from 1.34± l ;[gkeV to 0.44_^_keV. The 



-0.04 1 



electron plasma temperature kT e increased from 7 .6^S"| keV 
tol82_ 158 ' 



96 keV, whereas the optical depth of the plasma r p de- 
creased from 6.00!°'^ to 0.68_Jj_. This behavior reflects the 
known negative kT e and t p correlation (see Section 5.2). 



An Fe K n 



3.1. Fe line component 
emission component 



was detected at 



6.414+Q oof keV in the outburst data (see Fig. __ Inset), 
as well as at 6.42_q jj} keV in the 2nd observation (see Fig. [3] 
Inset). For the best fit values the centroid energies of XIS03 
and XIS 1 were left independent, due to a known shift in the 
energy calibration. Although the centroid energies for XIS 1 
generally show a slightly lower (30 e V) energy than the other 
XIS instruments^, our XIS 1 data instrument actually showed 
a 20-30 eV higher centroid energy. Residuals at ~ 7keV 
were observed and were modeled with an additional Gaussian 
line with a best-fit energy of 7.13 ±0.03keV, indicating the 
existence of an Fe line. The width a of the Fe line was 

15 http://www.astro.isas.jaxa.jp/suzaku/process/caveats/caveats_xrtxis06.html 
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set to the Fe K Q line value, and the normalization was left in- 
dependent. The obtained value for the normalization was 
consistent with the expected 12% of the K Q normalization. In 
the second observation the second line was very weak and its 
parameters could not be properly constrained. In this case the 
line centroid was frozen at the outburst value of 7 . 1 3 keV. The 
width was again coupled with the Fe K Q value. The residuals 
decreased from \ 2 = 721 for 477 dofs to the overall best 
fit values of 678 for 474 dofs, including two CRSFs, when 
the Fe line was included and an F-test showed only a 
2.8 x 10~ 6 probability that this improvement is of a statistical 
natureFl Another possibility for modeling the observed 
residuals at ~ 7 keV is the addition of an Fe K-edge, located 
at 7.1 keV. When including the edge, the improvement is only 
marginal (x 2 of 712 for 474 dofs) and the best fit energy is 
~ 7.6 keV. Freezing the energy to the value of 7.1 keV did 
not improve the fit. The remaining residuals indicated that an 
additional Fe K^j line would still be necessary. The measured 
energies of the Fe K Q and Kg lines were slightly higher than 
the expected laboratory values, but still in the same order of 
magnitude of the energy scale uncertainties of the instrument 
( ~ ± 20 eV @ 6keV according to the ABC guide)0 

Comparing the equivalent widths (EQWs) of the Fe K Q 
and K^j lines in the two observations showed no significant 
change during the outburst. The Fe K Q line showed an EQW 
of 51 ±4eV for both models, whereas the Fe line showed 
an EQW of 8 ± 3 eV for the cut o f f p 1 model and 10 ± 3 eV 
in the compTT model. For the second observation, the Fe K Q 
EQWs were 46 ± lOeV and 49 ± lOeV, and the K p EQWs 
were 8 ±2 eV and 11 ±9 eV for the cutoff pi and compTT 
models, respectively. The similarity of the Fe line EQW be- 
tween both observations indicates that the source of the Fe 
lines is near to the source of the continuum (see discussion). 

3.2. Cyclotron resonance scattering feature 

In the burst observation a strong residual at ~ 50 keV was 
observed in the HXD (Fig.U b). When including a CRSF-like 
feature with a Gaussian optical depth profile (gabs) the best 
fit was obtained with a centroid energy of 58 .2*2 - J keV, signif- 
icantly improving the x 2 from 1728 with 480 dofs to 752 for 
477 dofs for the cutof fpl model , confirming the di s cover y 
of the CRSF in the RXTE data by lDoroshenko et all (l2010t) . 
The width of the line was 14.1^ jkeV and the optical depth 
t was 60.1*j 5. Using the compTT model, similar values for 
£ C yc and a were obtained for the outburst data. In this case the 
optical depth was found to be 23.6t^% 4 - 

In the cutof fpl model, a second absorption line, using 
the gabs model, in the GSO data of the outburst spectrum at 
~ 1 lOkeV improved the x 2 to 678 with 474 dofs (see Fig. [2};). 
An F-test probability of 0.144 indicates that the line may not 
be significantly detected. However, an increase of the GSO 
Non X-ray Background (NXB) by 1 or 2% decreased the GSO 
counts to a non-detectable value. Another possible system- 
atic effect is the decay of the 153 Gd instrumental line, which 
is activated when the satellite passes through the south at- 
lantic anomaly (SAA), and thus creating a background line at 

16 See however Protassov et al. (2002) about the usage of the F-test in 
line-like features. 

17 An inflight energy calibration using the Mn-calibration sources of the 
XIS detectors could not be performed, due to the usage of the 1/4 window 
mode, which excludes the corner regions where the calibration sources are 
located. 



~ 150keV. This could impose a deviation at ~ 1 lOkeV if not 
accurately represented in the modeled background. Therefore 
caution is advised when interpreting this feature. 

In the lower luminosity observation only an indication of 
the fundamental CRSF was observed (see Fig. |3J. The line 
was fitted in 2 steps, where first the centroid energy and width 
were fixed to the outburst values, and the optical depth was 
left as a free parameter. With the addition of this CRSF, the 
cutof fpl model improved slightly from a x 2 /dofs of 605 
/ 432 to 574 / 431 with an optical depth of r = 15.4^;g. In 
the second step, the CRSF parameters were left free and the 
centroid energy decreased to 47 .4;^ 3 keV, whereas the opti- 
cal depth decreased to r = 6.0*2 1 an d tne width decreased 
to a = 5.Tf% keV. With a x 2 /dofs of 567 / 428 the best fit 
improved further (see Tab. [T]for details). To test the signif- 
icance of the CRSF in the second observation, 100000 PIN 
spectra were Monte Carlo simulated using the null hypothesis 
approach where a spectrum, including Poisson noise, was cre- 
ated using the best fit parameters without an CRSF line. The 
simulated spectra were then fitted with a continuum model 
and an additional gabs absorption line to test how often such 
a line could be detected in the spectral noise. The width was 
constrained to a value of a between 3.5 keV and 8.5 keV, so 
that neither very small features nor broad parts of the con- 
tinuum were modeled. Out of the 100000 simulations, 42% 
of the best fit results showed a non-zero value of r. When 
comparing the simulated r distribution with the best fit result 
of r = 6.0*1' j, ~ 96% of all simulated fits showed a value be- 
tween and the measured value. When taking the errors of the 
measured data into account (t < 3.9) the number is reduced 
to - 90%. 

For the compTT model, the addition of the fundamental 
CRSF with the frozen energy and width created residuals in 
the 10-20keV range, which could be significantly reduced by 
decoupling the optical depth of the compTT model in the PIN 
data from the XIS data. The best fit optical depth of the CRSF 
line was 27. 1+Jq changing the x 2 I dof values from 639 / 429 
to 567 / 428, respectively. When leaving the CRSF parame- 
ters free, the width of the CRSF increased to > 30keV, mak- 
ing it part of the continuum and tampering with the results. 
Freezing the CRSF parameters to the values obtained with the 
cutof fpl model did not result in a satisfactory overall fit. 

An alternative approach to describe the CRSF is to re- 
place the gabs m odel with the XSPEC model eye labs 
(Miha ra et al.| [l990). described by the resonance energy E Tes , 
the resonance width <7 les and the resonance depth r res . Using 
this model does not improve the fits significantly, and actu- 
ally results in a slightly worse x 2 /dofs value of 782 / 475 for 
the outburst data. Note that in the case of cyclabs the ra- 
tio between the energy of the fundamental and first harmonic 
line are fixed to 2, resulting in 1 more degree of freedom. 
The observed E KS of 49.5^ J keV is of the order of - 20% 
below the centroid energy obtained with the gabs model. 
This discrepancy stems from the use of a different calcula- 
tion of the line energy and is of the ord er of 10-20% lower 
than the measured gabs energy (see Nak aiima et al.l l20Toi 
for details). No significant changes in the continuum param- 
eters have been observed when using the cyclabs model. 
A width of cr ies = 19.4 +2 JkeV and a depth of r res = 1.45^;Jg 
are obtained for the outburst observation. In the second ob- 
servation the resonance energy shows a slight decrease (E KS = 
45.7*2 gkeV), and only slightly smaller than the gabs val- 
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FIG. 4. — Instrumental residuals at lower energies for XIS 1 (blue) and XIS 
03 (red) with the cutoff pi model, a) shows the clearly apparent lines are 
the Au K a line at ~ 2.1 keV and the Si K a line at ~ 1.8 keV. Additionally 
a possible Ni K Q line at ~ 0.9 which improves the fit slightly, b) shows the 
best fit residuals when the lines are included. 

ues. These results reflect again the difficulty of constraining 
the possible CRSF parameters in the second observation. The 
width and depth, cr res = 5.6+2 4 keV an( ^ r res = 0.46+j} are 
consistent with the gabs result showing a decrease in width 
and depth for both observations. 

3.3. Low energy calibration issues 

Strong residuals at lower energies, i.e. between 1.5 and 2.5 
keV (see Fig. were observed in both data sets. A compari- 
son with known background prop erties of the Suzaku/XIS in- 
strument (Yamaguchi et al. 2006) showed that these features 
are identical to the known instrumental Si K a and Au K Q 
lines, located at 1.74keV and 2.12keV, respectiv ely. In many 
obser vations of bright sources, e.g. 4U 19 07+07 dRivers et alj 
1201 01) and LMC X-3 dKubota et alj|2010l) , these energy bands 
are explicitly excluded. In this case, modeling these lines with 
two Gaussian emission features was sufficient to minimize the 
residuals and improve the x 2 /dofs from 1295 / 483 (without 
the lines) to the best fit value of 678 / 474. Note that the Au 
K Q line at ~ 2.1 keV could be described by the same Gaus- 
sian emission line at 2.21 ± 0.01 keV for both FI and BI XIS 
instruments, whereas the Si K Q line appears as an emission 
line for the BI XIS1 at 1.82 ± 0.01 keV (Fig.H blue) and a 
negative Gaussian line at 1 .89 ± 0.01 keV for the FI XIS 03 
combination (Fig. @] red), indicating that for the FI instru- 
ment the line is either over-subtracted or not properly energy 
calibrated, leading to a dip in the spectrum. By using two in- 
dependent lines at 1.89 keV (FI) and 1.82keV (BI) the residu- 
als can be well described. An additional Gaussian component 
at 0.93^;|] 2 keV, located close to the Ni K edge at 0.897 keV, 
slightly improves the residuals (x 2 /dofs of 731 / 477 with- 
out the line, compared to the best-fit value of 678 / 474) This 
improvement indicates a possible systematic error of the cal- 
ibration at lower energies for bright sources. This line is not 
very pronounced and can be omitted when using the compTT 
model. 

4. PHASE RESOLVED ANALYSIS 

For a phase resolved analysis the XIS and PIN data of both 
observations were folded with t he RX TE determined pulse 
period from iDoroshenko et all (12010b of P spin = 407.719 s, 

p spin = -4.6 x 10" 7 s/s and the MJD epoch of 54841.62. Pulse 
profiles with 20 phase bins were created for five different 
energy bands in the 1 - lOkeV energy range for XIS and 
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FIG. 5. — XIS and PIN pulse profiles for different energy bands for the 
outburst observation. PI, P2 and MEM indicate the regions used for the phase 
resolved spectral analysis. 

12-40keV energy range for PIN for both observations (Fig. [5] 
and Fig- [6j>. The statistical quality of the GSO data precluded 
the creation of pulse profiles and further spectral analysis. 

During the outburst, the pulse profile consisted of two 
peaks, where the main peak (PI, pulse phase 0.3-0.6) stayed 
dominant throughout all energy bands and became broader to- 
wards higher energies (Fig. |5). The second peak (P2, pulse 
phase 0.65 - 0.9) disappeared at energies above ~ lOkeV, 
consis tent with the RXTE observations by Doroshenk o et al.l 
(2010). The third region in Fig. [5] indicates the minimum 
of the pulse profile (MIN, pulse phase 0.95- 1.25), as deter- 
mined from the lower energies of the XIS instrument. 

In comparison, the pulse profile of the second observation 
(Fig. [6]) showed a similar Peak PI throughout the whole en- 
ergy band, although narrower at lower energies. No second 
peak was observed at the position of P2, although a small 
"bump" in the 5- 12keV energy range was still visible. At 
lower energies, i.e. for the 3-5 and 5- 12keV energy band, 
a small peak showed up on the opposite site of PI, the pulse 
phase where the minimum was defined in the outburst data. 
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FIG. 6. — Same as Fig.[3]for the second observation. The highlighted re- 
gions indicate the same phase bins as in the outburst observation for a direct 
comparison 

4.1. Phase resolved spectroscopy 

Spectra were extracted for three different pulse phases 
throughout the outburst: PI, P2 and MIN. XIS and PIN spec- 
tra were used and the same XIS grouping as in the phase av- 
eraged data was applied. The PIN data was again grouped to 
include at least 100 counts per spectral bin. The same spec- 
tral model was used as in the phase averaged analysis and the 
best-fit results are summarized in Table |2] In all pulse phases 
the CRSF component was visible in the spectra. Throughout 
the pulse phase the best-fit values for the galactic and intrinsic 
Nh values did not change significantly. Note that in both mod- 
els a slight decrease of the covering fraction during P2 can be 
observed. 

The most significant variation in the cutoff pi model 
was that of the cutoff energy E cut which decreased from 
29.5^ keV in PI, to 18.3^jkeV in P2 and down to 
10.7^'ij keV in MIN. At the same time the power law in- 
dex T varied from 0.42^ (PI), to 0.56^;$ (P2), and to 
-0.13^ for MIN. In the three phases the CRSF centroid 
energy changed from 66. 6^ '3 keV in PI, and declined to 



55.5^jkeV throughout P2, and 52.8!fgkeV for MIN. The 
CRSF width declined from 19.2^keV in PI to 11.4^|keV 
in P2, and to 9.3!^ keV in MIN. 

For the compTT model most parameters did not change 
throughout the pulse profile. The observed increase in the 
plasma temperature and the decrease of the optical depth r p is 
a known systematic anti-correlation (see discussion). As with 
the phase averaged data, the CRSF did not show the same 
changes throughout the pulse profile, as in the cutof fpl 
model. The best-fit values for the CRSF centroid energy were 
57.9!^ keV for PI, 57.7^ ( >keV for P2, and 51.0!^keV 
for MIN, and the widths were 13.5^ 5 4 5 keV,14.9!^keV, and 
9.0^ keV for PI, P2 and MIN, respectively. Fe K a and 
Fe energies were consistent with ~ 6.4keV and ^7.1 keV, 
similar to the cutof fpl values. 

5. DISCUSSION 

This paper presents an analysis of the two Suzaku observa- 
tions of the Be/X-Ray binary 1 A 1 118-61 during the peak of 
its outburst in 2009, January and ^13 days later. A CRSF, 
detected with RXTE at ~ 55 keV, could be observed in both 
observations, although the significance is lower in the sec- 
ond observation. An Fe line at 7.13keV has been ob- 
served in addition to the strong, narrow Fe K Q line at 6.4 keV. 
The broad band continuum was modeled with the empirical 
cutof fpl model, including an additional lOkeV systematic 
component to improve the residuals. Softening of the power 
law index T between peak and decay had been observed in 
an earlier outburst and could be confirmed. The Comptoniza- 
tion model compTT, where the lOkeV component was not 
needed, has also been applied. The pulse profiles at lower en- 
ergies changed from a two peaked to a single peaked profile 
between both observations. Phase resolved spectral analysis 
was performed for both observations and the same models as 
in the phase averaged analysis were applied. 

5.1. Outburst behavior 

The third observed outburst of 1A 1118-61 follows a pat- 
tern similar to the second outburst from 1992, i.e. a strong 
peak lasting ~ 3 weeks and an elevated level of emission up 
to 6 weeks afterwards. The time between outbursts was in 
both cases ~ 6200 days, corresponding to ~ 17 years, indi- 
cating that the outburst behavior could be periodic on very 
long time scal es. The prop osed orbital periods are ~ 35 days 
(ICorbetll 19861) . ~ 58 days (iReig. Fabregat & Cod 1 19971) and 
most recently 24 days established by Staubert et al.l (120101) us- 
ing the delay in pulse arrival time of RXTE monitoring obser- 
vations throughout this outburst. Using the latter method the 
Suzaku light curves are in full agreement with the RXTE re- 
sults (Staubert, priv. comm). A period of 24 days would put 
1A 11 18-61 in the wind accretor region on the "Corbet" dia- 
gram (ICorbetl l 1986). making it a very unique source for a Be 
system. 

A v ery similar scenario was introduced by Villad aet al.l 
(1999), who monitored the H a IHp emission before and af- 
ter the second outburst and proposed that the optical com- 
panion, Hen 3 - 640, has an extended large envelope where 
a weak interaction w ith the NS can occur. In the scenario of 
IVillada et"aT1 0999), the NS is orbiting the O star in an envi- 
ronment with gradually increasing density, until a steady ac- 
cretion disk is created and the X-ray flux suddenly increases. 
The sudden increase in the accretion material would provide 
the torque on the NS to produce the observed changes in the 
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TABLE 2 

Phase resolved outburst spectral parameters. Same models as in Table 1 



Parameter 






Cutoffpl 




Parameter 




CompTT 








Peak 1 


Peak 2 


Minimum 




Peak l 


Peak 2 


Minimum 


pnabs /VhLIU /cm A _ 




1 •> 5+0.02 

-0,03 


1.38 „„, 


1 71+004 
1 ' Z1 -0.04 




1 TQ+0 04 

1.28+»;»i 


1 m+U 6 


1 on+0 04 

12 16+™ 
-0,89 

70+ 04 
-Q.Q5 
ft S9+0 03 
u - -0.02 


pcfabs A^ H [10 22 /cm 2 


] 


11 l+° g 


n c:o+0 01 
0.58_ 002 


11.7+JI 




1 1 6+°* 


11.0+"' 


covering fract. 
blackbody kl [kevj 




66+™ 
-Q.Q3 
n ^9+0 o1 

ylD -0.02 


67+ 05 
-Q-A5 
n <:n+0.02 

°- 50 -0.02 




70+ 03 
n s/i+0 03 
°- 54 -0.02 


0.584? 

n ?n+0 03 
°- 59 -0.03 


blackbody norm [10 


-3 it 
J 1 


7 -°-1.6 


o a+0.6 
23 -0.4 


1+1.7 
S - 1 -0.7 




1 O -2 + 1.7 

12.3_ 19 


Q Q+1.6 

89 -1.6 


in 1+I.6 

10.1_ 14 


c cut [kevj 




TO c+7.2 

29.5,, 


IS 3+0.8 


in 7+1.3 


compTT To [keV] 


1 31+0.10 

1J -0.07 


1 41+0.04 


1 AS+0. 11 
L48 -0.07 

6.79+°^ 


r 




n j -,-1-0 06 


0.56«;W 


r, 1 -.4-0 OQ 


compTT kT [keV] 


7 n7 
'■-0.38 


11 /-\ 1 4-Q Ofi 

1 1 .91+2 77 






1 2+ 001 
u - lz -o.oi 


14+° 01 
u ' 1 ^-0.00 


o.o4-s : 88* 


compTT r p 


6 S1+ 31 
°- J1 -0.40 


4 00+° 47 

+ uu -0.45 


6.60!8; 3 f 










compTT norm t 


11+ 001 


05+° 02 
u - u;l -o.oi 


0.08+8»! 


Ecrsf [keV] 




66.6+°' 


55 5+0-5 






57.9+^ 


57.7+2.0 


51.0+4 9 


"crsf [keV] 
tcrsf 

EfeKa [keV] 




19.2+ 36 
130.8+g 6 


11.4^| 
44.1+2-3 
64 ,+"o 3 o 3 i 

°'^ 1 -0.Q! 


9.3+ 2 ' 5 
22 7+^6 

6 42+^ : ^ 
°-^ z -0.0l 




13.5+ l5:5 

64 ,+om 

°-+ 1 -Q.Ql 


149^:4 

79 j+A's 
ft 4 ,+(t OI 

°-^ 1 -Q.01 


9 0+4-6 

n s+5t?4 

1 D.J . 7 

°-+ z -0.01 


E Fc k/3 [keV] 




7 I6+ 005 
'■ lo -0.05 


7.08^ 






7 16+ 004 


7.o9+« : 8| 


7 1 7+0. 14 
'■ lz -0.12 


Eq. Width K a 1 Kg 


[eV] 


37+9/7+1 


59+1/11+9 


64+f/8!| 




39+1/10+' 


58+2/13+| 




Flux 2 -iokev[10 36 erf 


's/sec ]+ 


10.5 


8.14 


7.5 










CxiSl/Cpi N 




0.98/1.1 


0.97/1.04 


0.99/1.16 




1.04/1.23 


0.99/1.19 


1.08/1.23 


X 2 /dofs 




498/431 


508/420 


474 / 427 




528/438 


492 / 423 


542 / 430 



NOTE. — f Units in Photons keV cm " s , + unabsorbed flux using a distance of 5 kpc 



pulse period. The surrounding material is then swept out in 
a shor t time and the syste m returns to quiescence. Accord- 
ing to IVillada et all (fl999l) . an interval of 17 years between 
the outbursts is a reasonable time scale to accumulate enough 
material between the outbursts. 

5.2. Cyclotron features 

In the outburst, a CRSF has been observed at ~ 55keV 
and the parameters are consistent with the RXTE data 
dDoroshenko et a l. 2010). A CRSF in the second observation 
is expected but can not be confirmed with sufficient signifi- 
cance in the established data. 

CRSFs are generated by electrons with energies close to the 
ones determined by the discrete energy states of the Landau 
levels being excited to higher levels followed by de-excitation 
emitting a photon. In this process, photons from the accretion 
column above the magnetic poles with a resonant cyclotron 
line energy, are resonantly scattered out of the line of sight 
creating an absorption line-like feature. These features pro- 
vide a direct measure of the magnetic field strength close to 
the NS surface, where the fundamental cyclotron energy is 
given by 



p — 

'-'CMC — 



TieB 1 
m e c 1 +z 



1 1 .6 keV 

1+z 



xB 



12. 



(1) 



where B\2 is the magnetic field strength near the NS surface in 
units of 10 12 Gauss and z is the gravitational redshift. Assum- 
ing a typical NS mass of 1.4M Q and NS radius of 10 km gives 
z = 0.3. For 1 A 1 1 18-61 the CRSF at ~ 55 keV corresponds 
to a magnetic field of ~ 6.4 x 10 12 Gauss, making it, together 
with the ~ 50 keV line of A 0535+26, one of the strongest ob- 
served magnetic fields on an accreting neutron star in a binary 
system. 

A possible feature in the residuals at ~ HOkeV could in- 
dicate the existence of a first harmonic line, as observed in 
multiple sour ces, such as 4U 01 15+63 dHeindl et aljjl999b . 
4U 1907+09 (Cus umano et alj|1998t iMakishima et all 1 19991) 
and VelaX-1 (Makishi ma et alj Il999t iKre vkenboh m et al] 



2002). Including a gabs line at this energy does not sig- 
nificantly improve the fit, though a 153 Gd inst rumental line at 
~150keV, which is due to electron-capture (Kokubunet al. 
120071) . as well as possible systematic uncertainties in the back- 
ground weaken the significance of the detection further. 

Assuming that the observed width of the CRSF is due to 
the Doppler broadening of the electrons responsible for the 
resonances, one can estimate the inferred plasma tempera- 
ture in the CRSF region using equation 4. 1 .28 from M eszarosl 
(fl992h : 



Awd = WCyc 



2kT 



1/2 



|cos(9)|, 



(2) 



where Awd corresponds to the CRSF FWHM which is cal- 
culated from crcyc, wcyc is equivalent to the CRSF centroid 
energy £cyc an d @ is the angle between magnetic field lines 
and the line of sight. With cos(8) = 1 the lower limits for the 
plasma temperature are kT k 20.8 keV for the cutoffpl 
model and kT « 12.6 keV for the compTT model when using 
the values from Tab. 1 . 

Although the existence of a CRSF in the second observa- 
tion is only marginal, it is very likely that such a line ex- 
ists at lower luminosities. Luminosity dependance of CRSFs 
have been observed in multiple other binary system. In Be- 
Binari es, such as V 033 2 +53, 4 U01 15+63, and X 03 3 1+53 
where Tsygan kov et al.l (12006b and [Nakajima e t al .1 (120061 
2010) found a negative correlation between the CRSF cen- 
troid energy and the luminosity of the s ource for luminosities 
approaching the Eddington luminosity. Stau bert et al.l (120071) 
on the other hand found a positive CRSF-luminosity corre- 
lation in the low mass X-ray binary HerX-1 for luminosities 
far below the Eddington luminosity. 

The type of correlation seems to depend on whether the ob- 
served luminosity is above or below the Eddington luminos- 
ity. When the luminosity is above the Eddington luminosity, 
the infalling protons start to interact before they are part of 
the accretion column and a "shock front" region, where the 
CRSF most likely occurs, is created. With increasing lumi- 
nosity, the proton interaction occurs farther away from the NS 
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surface, where the magnetic field is lower and therefore the 
observed CRSF is seen at lower energies. Below the Edding- 
ton luminosity, the proton interaction does not occur above 
the accretion column but is part of it. When the luminosity 
increases, the accretion "pressure" increases as well and the 
region where the CRSF is created is pressed closer to the NS 
surface, where the magnetic field is higher. This results in 
a positive CRSF-luminosity correlation, as observed in Her 
X-l. 

The flux levels (Tab. 1) obtained for both observations 
showed that the first observation is most likely above, and 
the second observation is definitely below the Eddington lu- 
minosity of the system, so that an anti-correlation between 
CRSF centroid energy and luminosity are expected. 

5.3. Continuum comparison 

When comparing both observations, changes in the broad 

band continuu m parameters were observe d. During the 

1974 outburst lives, Sanford & Bell Burnetii (119751) observed 
a harder spectrum in the peak compared to later observation. 
Looking at the cutoff pi model, the Suzaku data showed 
a similar behavior with a very hard spectrum (power law in- 
dex r = 0.35) at the luminosity peak, and a much softer power 
law index of T ~ 1 in the second observation. Note also that 
the cutoff/folding energy increased with declining luminosity. 
A similar behavior has been observed in a number of differ- 
ent HMXB trans ients, such as A 0535 +26 dCaballerol 120091) 
and V 0332+53 dMowlavi et all 120061). This is in contrast 
to EXO 2030+375 dReynolds.' Parmar & Whitelfl993l). where 
the observed folding en ergy increased with lower luminos- 
ity. Upong et al] (119901) interpreted this parameter in phase 
resolved results and concluded that the folding energy reflects 
a change of the viewing angle on the accretion column, which 
allows a deeper look into the emission region and therefore 
directly correlates with the observed electron plasma temper- 
ature. The change observed in the phase averaged observa- 
tions follows a similar reasoning, and at lower luminosities 
the emission region is closer to the NS, where the observed 
plasma temperatures are expected to be higher. 

A co mparison of comp TT parameters between Suzaku and 
RXTE (Dorosh enko et al.l2010h showed consistent results for 
the spectral parameters in the outburst. In contrast to the 
RXTE results, the observed column density Nu is higher and 
consistent with the results obtained by the cut o f f pi model. 
The main difference between these results and the RXTE data 
lies in the use of a combined Galactic and intrinsic column 
density and the additional need of a low energy black body, 
so that these values cannot be com pared directly. Note that 
the RXTE/PCA instrument used in iDoroshenko et all (120101) 
is not sensitive to data below 3keV and therefore a partial 
covering, as well as a black body component could not be 
modeled. Both, the black body temperature for the low energy 
excess and the photon seed temperature in the compTT model 
decreased by ~ 1/3 between outburst and the second obser- 
vation, indicating a correlation between both parameters. A 
possible explanation is that the source of the compTT seed 
photons is the same as the soft excess modeled by the black 
body component. The observed change in the optical thick- 
ness indicates the the observed material is optically thick in 
the outburst and gets optically thinner for the second observa- 
tion. 

Although an increase in the plasma temperature of the 
compTT model was observed after the outburst, i.e., mainly 
consistent with the interpretation of the cutoff pi param- 



eters, one has to be careful with a direct interpretation of 
the spectral parameter. The optical depth r and the electron 
plasma t emperature kT R sho w a very strong negative correla- 
tion (e.g. Wilms et al. 2006), and the best-fit values cannot be 
used directly for interpre tation. Following the definition in 
iRybicki & Lightrrtanl (119791) . the Compton y parameter: 

4-kT 

V= |max(r,r 2 ), (3) 

can he lp in the physical interpretation. iReynolds & Nowakl 
d2003l) used the Compton y parameter in the description of 
accretion disk coronae around black holes. A value of y 1 
or slightly higher means that the average emitted photon en- 
ergy increases by an "amplification factor", A(y) m exp(y), and 
is referred to as "unsaturated inverse Comptonization". For 
y 1 the average photon energy reaches the thermal energy 
of the electrons. This case is called the "saturated inverse 
Comptonization". Using kT and t p results in y = 2.13 ±0.2 
for the outburst data and y = Uq 7 for the second observa- 
tion. Typical calculated values of y-parameters are smalle r 
than 1, e.g., ~ 0.5 for CvgX-1 dRevnolds & Nowakl l2003l). 
~ 0.2 for 4U 2 206+54 dTorreion et al.l 12004b. and~^ 0.6 for 
XV 1832-330 dParmaretal.1120011) . Note that the second ob- 
servation is very badly constrained due to the large error bars 
on the electron plasma temperature and the optical depth and 
therefore cannot be used for interpretation. The y parameter 
for the Suzaku outburst indicates that the system is in, or very 
close, to a "saturated inverse Comptonization" state. 

5.4. Fe lines 

In both observations an Fe K Q and Fe emission line has 
been observed, where the FeK Q - FeK^ normalization ratio 
of 12% is consistent with neutral material. The FeK Q EQWs 
were 51 ±4eV and 46 ± 10 eV for the first and second ob- 
servation, respectively, while at the same time the 5-7 keV 
power law flux dropped from 6.50 ±0.01 x 10~ 10 ergs cm" 2 
s _1 to 1.23 ±0.01 ± xlO" 10 ergs cm -2 s _1 . The relatively 
constant EQWs imply that the Fe line emitting region is rel- 
atively close to the source of ionizing flux, so that the ob- 
served Fe line intensity adapts quickly to the changing inci- 
dent flux. The Fe line normalization was 1 .79 ± 0.06 x 10~ 3 
photons cm" " 2 s" 1 and 3.15 ±0.3 x 10" 4 photons cm s for 
the respective observations. The ratios between Fe normaliza- 
tion and power law flux were consistent in both observations, 
2. 75 ±0.10 and 2. 56 ±0.27 x 10 6 photons/ergs, reflecting the 
two calculated EQWs. Using the continuum flux difference 
and the time between observations, the calculated change in 
power law flux was ~ 0.415 ±0.001 x 10~ 10 ergs cm" 2 s -1 per 
day in the 12.7 day period between the observations. A linear 
decrease could be assumed using the BAT data of Fig. 1. For 
the second observation, the Fe normalization was manually 
increased in XSPEC until the resulting EQW matched the up- 
per limit of the measured EQW. This value corresponds to an 
Fe normalization of ~ 3.9 x 10 -4 photons cm" 2 s , assuming 
the measured continuum given above for the second obser- 
vation. Together with the previously calculated constant of 
2.56 x 10 6 , a matching power law flux of ~ 1 .52 could be cal- 
culated. Using this flux difference, together with the rate of 
change in the 5-7 ke V energy range, one can then determine a 
maximum time delay which would still preserve the observed 
EQW within errors. The calculated upper limit to the delay 
between the Fe and X-ray emission region is ~ 0.7 days. Note 
that this estimation is very simple and uncertainties from the 
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TABLE 3 
Pulse fractions for 
1 A 1 1 1 8-6 1 



Energy [keV] Outburst 2nd Obs 



1-3 


0.26 


0.73 


3-5 


0.31 


0.90 


5-12 


0.48 


0.53 


12-20 


1.00 


0.78 


20-40 


1.19 


0.83 



orbital motion or from line of sight assumptions are not taken 
into account here. 

5.5. Phase resolved description 

1A 11 18-61 shows a similar energy and luminosity pulse 
profile dependency as m any other so urces, e.g. 4U 01 15+63 
ilTsygankov e"tal]l2007h. V 0332 +53 dTsygankov et al.ll2006h 
and A 0535+26 dCaballerol2009l) . In the outburst observation, 
the observed broad double peaked structure at lower energies 
changed to a single peak profile above 10 keV, where the main 
peak (PI) broadened slightly with increased energy and the 
secondary peak (P2) weakened. For the second observation, 
the pulse profile changed significantly. There was still a dom- 
inant primary peak at PI, although it is narrower in the lowest 
energy band. The secondary peak is only marginally indicated 
between 5- 12keV (see Fig. [6) and absent in the other energy 
bands. On the other hand, between 3 - 12keV a small peak or 
extended shoulder could be observed preceding the primary 
peak. 

To confirm the disappearance of the second peak, pulse 
fractions of the total counts of the P1/P2 regions were cal- 
culated for the indicated energy bands. The pulse fraction 
was defined as (Pl C ountrate-P2 C ountrate) / P2 C ountrate where the PI 
and P2 regions are indicated in Fig. [5] and [6] Table [3] shows 
that the pulse fraction increases towards higher energies in the 
outburst data, typical for a weakening of the second pulse. In 
the second observation the values do not vary for the differ- 
ent energy bands, and the small peak in the 5- 12keV band 
shows only a marginal smaller pulse fraction than in the other 
e nergy bands. 

Tsygank ov et al.1 ((2007) described similar profiles and pro- 
posed that a misalignment of the rotational and magnetic axes 
of the neutron star leads to the case where one accretion col- 
umn is observed whole, whereas the second accretion column 
is partially screened by the NS surface so that only the softer 
photons, which are created in the higher regions of the ac- 
cretion column, are observed. When the overall luminosity, 
and therefore accretion rate, decreases, the column height de- 
creases and even the soft photons from the second pole are 
shielded. This behavior is in contrast to the variation observed 
in the cyclotron line production region, where the line form- 
ing region is closer to the NS surface at higher luminosities . 

For a more physical picture gravitational effects, such as 
light bending, h ave to be taken in t o account when discussing 
pulse profiles dKraus et al.1 1 19951: iMeszaros & Nagell 11985b . 
For a canonical NS with a mass of 1.4M© and a radius of 
10 km the visible surface of a NS is 83%. With this increased 
surface visibility the pulsed flux from both hot spots is visible 
over a longer part of the pulse phase, and parts of the accre- 
tion column of the second peak are still visible, although that 
hot spot is on the far side of the NS. With decreasing luminos- 
ity, the hot spot size decreases and the visible fraction of the 
second accretion column disappears. 



Pulse profile decompos ition methods have been devel- 
oped by iKraus et all d 19951) . and ha ve been applied to mul- 
tiple source s e.g. EXO 203 0+375 dSasaki et al.1 120101) and 
A 0535+26 dCaballeroir2011l) . Under the assumption of a 
slightly distorted magnetic dipole, with this method it is pos- 
sible to disentangle the contribution of the two emission re- 
gions, and constraints on the geometry of the pulsar and on its 
beam pattern can be obtained. The application of this method 
on the 1 A 1 1 1 8 data presented here will be left for future work . 

A more physical a pproach for a b r oad ba nd spectrum has 
been introduced by Becker & Wolff (2007), using bulk and 
dynamical Comptonization of photons in the accretion col- 
umn. First results on 4U 01 15+63 have been promising 
dFerrigno et al.l2009l) and 1A 1 118-61 is a good candidate for 
future tests. ISchonherr et al1 d2007) developed a new physical 
model for the CRSFs based on Monte Carlo simulations for 
the Green's functions for the radiative transport through a ho- 
mogeneous plasma, which had b een successfully applied for 
different sources, e.g., CenX -3 (ISuchvet al.l 12008b . Due to 
advances in the overall code dSchwarmll2007h the test of this 
model on this data set is beyond the scope of this paper. 

6. SUMMARY AND CONCLUSIONS 

In this paper we analyzed the broad band spectrum and 
pulse profiles of 1A 1118-61 during the peak of its third 
observed outburst and compared the results with a second 
observation which occurred ~ 2 weeks after the main peak. 
The time between outbursts is consistent with a continuous 
low le vel accretion mechanism as suggested by Villa da et al.l 
(1999). A CRSF during the outburst has been confirmed at 
~ 55 keV indicating one of the highest known B-fields ob- 
served in HMXBs. In the second observation there is only 
a weak indication of a CRSF fundamental line. A change in 
the CRSF centroid energy with respect to luminosity would 
be expected and would help to understand the physical en- 
vironment close to the NS surface. Variations of the CRSF 
can also be observed in the phase resolved spectroscopy of 
the outburst observation. The calculated y parameters show 
that the inverse Comptonization during the outburst is very 
close to saturation. The ratio between the FeK Q and FeK,3 
normalization of ~ 12% shows that the emitting material is 
mostly neutral. Using the compTT model, one can deduce 
that the emitting material is optically thick during the out- 
burst and optically thin in the second observation, although 
the known kT — r correlation has to be considered. The pulse 
profile change for different energy bands and luminosities is 
similar to other observed HMXBs and can be explained with 
a misalignment of the rotation axis to the magnetic field axis. 
A change in the pulse profile shape with luminosity has been 
observed and indicates that the visibility of the second hot 
spot changes between the observations. Phase resolved anal- 
ysis throughout the outburst also indicates a change in the 
observed magnetic field, which could be caused by differ- 
ent viewing angles onto the accretion column. Future work 
with a pulse deconvolution technique, as well as a more phys- 
ical model will provide a better understanding of the involved 
physical processes. 
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